Final State Interactions in Decays of the Exotic tti Meson 
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We analyze the role of final state interactions (FSI) in decay of the lighest exotic meson, with 
J^^ — 1 ^ TTi. We use the relativistic Lippmann-Schwinger equation for two coupled Trfei and 7rp 
channels. The first one is the predicted dominant decay mode of the tti, whereas in the other a 
narrow 7ri(1600) exotic signal has been reported by the E852 collaboration. The FSI potential is 
constructed, based on the lj meson exchange between the two channels. We find that this process 
introduces corrections to the tti widths of the order of only a few MeV. Therefore, we conclude that 
a substantial irp mode cannot be generated through level mixing. 
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I. INTRODUCTION 

Exotic mesons may provide a unique experimental 
handle on gluonic excitations. By definition, spin ex- 
otic mesons have quantum numbers, spin (J), parity 
[P) and charge conjugation (C), which cannot be ob- 
tained by combining quantum numbers of the valence 
quark and antiquark alone. Thus other degrees of free- 
dom, in addition to the valence components, must be 
present in the leading Fock space component of an ex- 
otic meson to generate exotic combinations of J^'^ , e.g. 
J^*-^ = 0~~,0"' ,1 '",2 Since multiquark states 

are expected to fall apart, it is expected that the lead- 
ing components of exotic mesons will contain excitation 
of the gluonic field in the presence of just the valence 
QQ pair. This is confirmed by lattice simulations [H Q 
which find that the lightest exotic meson has J^'~^ = 1 ^ 
and overlaps with the state created from the vacuum by 
application of operators with Qj Q E^nd gluon field (link) 
operators. The mass of this lightest exotic is expected be- 
tween 1.6 — 2 GeV 0,^1^1^] with the uncertainty in the 
lattice estimates coming from chiral extrapolations 0- 
There are also estimates of exotic meson properties based 
on the Born-Oppenheimer approximation which indicate 
that the gluon, if treated as a constituent particle in the 
ground state exotic meson, has to carry one unit of or- 
bital angular momentum with respect to the QQ 8j . 

Whether exotic mesons have been seen experimentally 
is an open question. Some tantalizing candidates exist 
in particular in the ti'tt and pir decay channels 9, 1^ . 
The J^'-' = 1 ^ partial wave in the 77' tt system mea- 
sured in TT^p rj'TT~p is as strong as the J^'-' — 2+"'" 
wave from the 02 resonance decay. Unfortunately, the 
exotic wave peaks at M^/^r '--^ 1.6 GeV where there are no 
other well established resonances decaying to rj'ir. There- 
fore extraction of the phase of the exotic wave, needed to 
establish its dynamical character, depends on model as- 
sumptions regarding the behavior of other waves. With 
the GlueX/Hall D experiment one should be able to cir- 



channel study together with the analysis of production 
mechanisms. If current experimental results do indeed 
correspond to elementary QCD exotic mesons, then there 
might be a potential discrepancy between model pre- 
dictions and the experimental data. Most models pre- 
dict that exotic mesons should primarily decay to "ex- 
otic" channels dominated by two meson final states with 
one of them being an S-wave and the other P-wave me- 
son HI [H [H [m [H e.g. 5i7r, /itt, Ki,it, while the 
experimental sightings come for analysis of " normal" de- 
cay channels i.e. with two S'-wave mesons, like the 77' tt 
and the p-n modes. In this paper we examine the effects 
of final state rescattering between mesons as a potential 
mechanism for shifting strength from one channel to an- 
other. 

In a previous paper |l5j | we constructed normal and 
hybrid meson states and studied kinematical relativistic 
effects at the quark and gluon level. In this work we will 
estimate the size of corrections to the tti decays origi- 
nating from the meson exchange forces between mesons. 
Since particle number is not conserved and the parti- 
cle momenta are of the same order as their masses, this 
problem should be treated in a relativistic formalism. In 
Sec. II we apply the Lippmann-Schwinger equation for 
the two coupled tt/j and 7r6i channels. The first is the 
channel in which a narrow 7ri(1600) exotic has been ob- 
served H3|; ^iid the latter is the predicted dominant 
decay mode of the tti. Then in Sec. Ill we describe the 
computational procedure of solving the resulting integral 
equations and we present a discussion of the numerical 
results. 



II. FINAL STATE INTERACTIONS IN tti 
DECAY 

There exist several models of exotic meson decays, and 
all of them predict a small branching ratio to the 7rp 
channel pdl Ha. Il4| . It is possible, however, that this is 



2 



where 



FIG. 1; Final state interaction nbi ^ up. 



G^{E) - [E-Ho{i)+ie] 



(4) 
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FIG. 3: Diagrammatical representation of the Lippmann- 
Schwinger equation. A bold vertical solid line corresponds 
to the total amplitude of the interaction between two-meson 
states, whereas a normal vertical solid line represents the sum 
over all intermediate states. 



FIG. 2: The lowest order FSI corrections to the tti decay 
from the interaction vrfoi ^ up. A bold horizontal solid line 
represents a hybrid and normal horizontal solid lines refer 
to mesons. A vertical dashed line corresponds to a single 
uj meson exchange. The circle represents the overall decay 
amplitude. 

can subsequently decay into tt and lo, and the lo can ab- 
sorb the other tt and produce a p, as shown in Fig.^ The 
w exchange is expected to dominate the mixing potential 
since this is the dominant decay mode of the 6i meson. A 
similar mixing effect (mediated by p exchange) in connec- 
tion with the rjiT P-wave enhancement in the 1400 MeV 
mass region was considered by Donnachie and Page [l6l | 
in perturbation theory. However, since the meson cou- 
plings are large, one must sum up all possible amplitudes; 
with the lowest order diagrams shown in Fig. [21 In order 
to describe the total contribution of the final state inter- 
actions to the decay width of the tti , we need to solve the 
Lippmann-Schwinger equation, 

T = V + VGT, (1) 

which sums up the ladder of w exchanges and the bare 
exotic meson exchanges. 

We will denote the single particle exotic state tti by 
I a), and the two-particle states irp and tt&i by Roman 
letters. We first introduce the matrix elements for the 
potential V, 

V^, = = (2) 

and similarly for T. The elements Vai are just the am- 
plitudes of the corresponding " bare" decays of the tti de- 
scribing transitions at the quark level, whereas Vij are re- 
lated to the final state interaction potential. In our state 
space the Lippmann-Schwinger equation can be written 
as. 



+ I I + I I I + - 



FIG. 4: The Lippmann-Schwinger equation in a meson sub- 
space. 

Both equations in Eq. are represented diagrammat- 
ically in Fig.|21 Hereinafter, if an index appears twice or 
more then a summation over this index is implicitly as- 
sumed. The Hamiltonian of the state \i) is Hf){i), e — > 0+, 
and E is the energy which will be taken as the mass of 
the TTi . The matrix elements can be eliminated using 
the second equation in Eq. (|3J), leading to 

Ta,^ - Vc,, + V^,G, [1 - VG]-T^ {Vk^ + Vk^Go^T^i). (5) 

Now we introduce a subset of the T matrix which acts 
only between states \'Kp) and |7r6i). This is denoted by t 
and is defined by 

t^V + VGt, = [1 - VG]r^' Vkj . (6) 

Diagrammatically this represents the sum of all diagrams 
without hybrid intermediate states, as shown in Fig. ^ 
Therefore we obtain 

Taj — ^cti ^otjGjtji -\- 

+ VajGj[Sjk + {tG)jk]VkaGaTai — 

= iV[l - GV]-%, + [Vil + Gt)GV]^^G^T^,. (7) 

Solving for T^i gives, 

Ta^ = G^MG^' - Sj-i(y[l - GVr%,, (8) 

where Eq = {V{1 + Gt)GV)aa is the self-energy of the 
TTi, shown in Fig. [S] Both and Sq, are c- numbers, 
and thus 
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E = 




FIG. 5: The hybrid self-energy. 




FIG. 6: FSI correction to the vri decay amphtude. 

where we again used [1 — GF]~^ = 1 + Gt. A resonance 
production process has the form of z ^ ttiX — > M1M2X, 
where i denotes an initial state and X denotes possible 
other particles in the final state, assumed not to interact 
with the resonance tti, and Mi (M2) is the final state 
meson (isobar) state. The full amplitude has the form 

A{i TTiX M1M2X) = PaGcTct = 

= Pa-^ i —{V + VGt)^,, (10) 

I 



where Pa describes the coupling between the bare exotic 
and the initial state. Thus in the full amplitude the bare 
exotic propagator [E — raex)~^ is replaced by the dressed 
one, {E — rriex — ^q)~^- Furthermore if S„, which is 
energy dependent, is expanded around the physical exotic 
mass, one obtains the Breit-Wigner propagator. Here, 
we are interested in the effective coupling of the exotic to 
the two-meson channels. This is determined by the last 
term in Eq. H1U|I . In the absence of FSI, the coupling is 
determined by the matrix V which changes toV + VGt 
due to re-scattering. We will therefore be studying the 
reduced T matrix element. 



with tji defined by . The diagrammatic representation 
of Eq. is given in Fig. |H1 

Eq. lO written in a full notation in the rest frame of 
the TTi is given by 



i(p,q,A,A') = V(p,q,A,A') + 5] / — -^^--T^-T-l/(P,k,A,A")G(fc)t(k,q,A",A'), (12) 

(27r)-^4tJi(fc)w2(fc) 

I 



where 

uji{k) = £;(TOi,k), Ho{k) =uji{k) +uJ2ik), 
G{k) ^ [E - Ho{k) + (13) 

and rrii, i — 1,2 are the masses of mesons in the inter- 



mediate two- meson state (related to G). In the above p 
and q are the relativistic relative three-momenta between 
two mesons, and the dependence on the center-of-mass 
momentum has been already factored out. Spins A refer 
to either bi or p. 

We introduce the partial wave potentials by 



VLL'{p,q)^ J2 dnpdn^{L,M;l,X\J,j){L',M'-,l,X'\J,j)V{p,(i,X,X')YLM{p)Yl,M,{q), (14) 

M,M',X,X',j •' 

I 

where dfJk is the element of the solid angle in the di- tLL'ip,q)- Substituting VllKPi?) i^to (|12|l gives, 
rection of the vector k and k = |k|. Similarly we define 



tLL' (p, q) = Vll' (p> 9) + X] / (27r 



k^dk 



(27r)34cui-.,i(fc)wi..,2(fc) 
I 



VLL'ip,k)GL"ik)tL'^L'{k,q), 



(15) 
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In our state space we can have L — 0,2 (the relative 
angular momentum between tt and 61) or L = 1 (between 
TT and p). For a tti we also have J = 1. Thus in Eq. H15|) 
we must substitute: 

mL.i^myr, m.0,2 = m-2,2 = TObi , mi,2=mp. (17) 

Parity conservation reduced the number of non- vanishing 
matrix elements of the final state interaction potential to 
Vol, Vio, V12 and V21. Moreover, from CP invariance we 
have Vqi = V{q and 1^12 = 1^2*1 ■ The integral equation l(T^ 
cannot be solved analytically and one needs to replace it 
by a set of matrix equations. The details will be given in 
the next section. When this is done and III' are found, 
we may go back to Eq. Hll|l which becomes 

v-^ f k^dk 
~a,iP) --dP) + Y.j (2.)3(2J+l)4.,(^V2(fc) 

xaL'{k)G{k)LdL'L{k,P), (18) 

where aL^P) are the partial decay amplitudes defined in 
Ref. I'^ order to obtain the corrected widths, they 

must be replaced by the corrected amplitudes aL{P). 

Finally, we proceed to the form of the final state in- 
teraction potential. The dominant effective interactions 
between the &i7r and the pir channels are expected to 
originate from u) exchange. The 61 — > wtt is the domi- 
nant decay channel of the hi meson and the pun: coupling 
is also known to be large. The effective Lagrangian for 
the p-KLO vertex is 

Lp^^ = ffp.c.e''"^"a^w.7r^aAp;, (19) 

and for the bi'KLj vertex is 

Here the index i corresponds to isospin. The coupling 
constant, gp^,^ lies between 0.01 MeV"^ and 0.02 MeV^^ 
and we will take a value 0.014 MeV~^ 0, whereas gbnruj, 
can be obtained by calculating the width of the decay 
bi TTLo and comparing with its experimental value of 
142 MeV 18|. Consequently, the amplitude for 61 ttu 
is equal to 

A = -gb,^^e\Xb,)e"{X^,P), (21) 



where P is given by 

E{m^,F)+E{m^,P)=mb„ (22) 

and e*(A) are the spin-1 polarization four- vectors. The 
absolute value of the square of this amplitude, summed 
over and averaged with respect to Xb^ , gives 

\A\' = glUl + ^l (23) 

and using Ts-wave = 142 McV- 0.92 = 131 MeV (the fac- 
tor 0.92 comes from the partial wave D/S ratio) we get 
9bi-^uj = 3650 MeV. Therefore gpsi = 9p-ruj ■ 5bi7rt^ is ap- 
proximately 50. 

The final state interaction potential can be obtained 
from Lagrangians of Eqs. H19|l and H20|) expressed in 
momentum space and dressed with the instantaneous lo 
propagator, 

V{p, q, Xb, , Xp) = gpSI^t^yarP^q"-, 4" 2 

(p - q)^ + ml 

xe-(A,,,-p)e"(Ap,q), (24) 

where p is the momentum of the vr in the |7rfei) state and 
q is the momentum of the p. At short distances, (for 
large values of p and q) this potential is singular which is 
a consequence of treating mesons as elementary particles. 
Therefore we must regulate this potential with an extra 
factor that tends to zero for large momenta. We will 
choose an exponential function 

e-IP-l/^, (25) 

where A is a scale parameter expected to be on the order 
of the inverse of the meson radius A ~ 0.5 — 1 GeV. 



III. NUMERICAL RESULTS 

The Lippmann-Schwinger integral equation of Eq. I|15|) 
corresponds to outgoing wave boundary conditions. This 
means that the singularity of the term G{k) is handled by 
giving the energy E a small positive imaginary part ie. 
An integral of this form may be solved using the Cauchy 
principal-value prescription. 



f{k)dk 



E-Ho{k) 



f{k)dk 



i7r/(fco 



fdHoik) 



dk 



k^ 



1,2 



/(fc)(fc2 



k^ 



E-Ho{k) ^■■•'^■-'V 
I 2f[ko)LOi{ka)L02{kn 



E - Hoik) 



E 



dk 

- i7r/(fco 



k—kQ 

^l(fco)^2(fco) 
knE 



(26) 



r 
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gularity: 



tLL' {p, q) = Vll' {p, 9) + X! 



dk 



i{2ny{P-kl^ 



in 



E 



[S-i/o,L"(fc)KM(fcVL",2(fc) 

4(27r)3£; 



(28) 



where 

Ho,L{k) = WL,i(fc) + ujlA^) 
and the quantities fcp.L are defined by 
Hox{ko,L)--E. 



(29) 



(30) 



In the decay tti nbi, the D-wave amphtude is neg- 
hgible compared to the S-wave, thus we may reduce our 
angular momentum space to L = 0, 1. Therefore, we 
only need the formula for Voi{p, q), and after expressing 
Clebsch-Gordan coefficients in Eq. (|14|) in terms of the 
spin-1 polarization vectors we obtain 



K)i(p, q) 



in 



I df^prfl7q^(p,q,A,A') 

A, A' 



(31) 



where T^(p, q, A, A') was given in Eq. (|24|l . The S-matrix 
is obtained from the t matrix via 



.y/ko,Lko,L' , , 
i>LL' = OlL' - « Tr^TF. tLL'{kQ,L,ko,L')- 



(32) 



As a 2 X 2 unitary matrix it can be parametrized by two 
real scattering phase shifts So and Si , and one inelasticity 
parameter, ij, (0 < 77 < 1), 



Sll' — 



2iSo 



jy2gi(5Q+i5i) 



(33) 



Eq. may be solved, as we already mentioned, by 
converting the integration over k into a sum over N in- 
tegration points fc„, n=l,2,...A^ (determined by Gaus- 
sian quadrature) with weights Wn |2(]|. 

Let us assume that the coupling constant gpsi is a 
variable, and we introduce the potential strength / = 

gpsi/dpsi^ where gp^j = 50. In Figs.[7|and|Hlwe present 
the phase shifts as functions of /, for the energy E = 1.6 
GeV and the scale parameter A = 2 GeV. Their depen- 
dence on i5 for / = 1 and A = 2 GeV is shown in Fig. |51 
Within an experimental range of the expected hybrid 
mass rriex — E and the hadrodynamical coupling con- 




FIG. 7: Phase shift So in degrees for the interaction 7r6i — Tip 
as a function of the potential strength I, for E — 1.6 GeV 
and A = 2 GeV. 




FIG. 8: Phase shift Si in degrees for the interaction vrbi — -Kp 
as a function of the potential strength J, for E — 1.6 GeV 
and A = 2 GeV. 



niex is large, thus the Born approximation is accurate. 
Furthermore, smaller values of closer to the nbi 

threshold, which (together with the second power of the 
FSI potential in momenta) should make the phase shifts 
small. Therefore, one expects that the final state in- 
teraction should not dramatically change the widths of 
TTi —t nbi and tti np obtained without final state in- 
teractions, such as was discussed in Ref. [Tsf. 

In Tables and |n] we compare the original widths (in 
MeV) with the FSI-corrected ones for various values of 
m^x and A, for 1 = 2. The strength was doubled in order 
to compensate for the effect of the regulating factor in 
the FSI potential. For the values of A other than 2 GeV, 
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separable potential, 



(34) 



14UU 161X1 I7UU ISIHI 1900 21)U0 

ra JMeV 

FIG. 9: Phase shifts So and 5i in degrees, for the interaction 
vrfoi — TTp as functions of the energy E, for I — 1 and A = 2 
GeV. 

qo are the relative momenta between mesons for the nbi 
and TTp channels, respectively. In Figs.llOlandllllwe show 
their dependence on / for niex = 1-6 GeV and A = 2 
GeV. The behavior of 6i indicates the existence of np 
resonances for the free parameters given. 

From Figs. ^1 and ^2 it follows that the irbi chan- 
nel can change significantly, if we increase the potential 
strength enough. However, even for large potentials the 
TTp channel width remains on the order of a few MeV. 
These results, which are shown in Tables and |nl are 
only weakly dependent on the value of the parameter A. 
Therefore, final state interactions cannot produce a large 
width for the tti Tip mode. 



niex [GeV] 


no FSI 


A = 0.5 


A = 1 


A = 2 [GeV] 


1.4 


85 


84 


86 


87 


1.5 


153 


149 


151 


153 


1.6 


150 


144 


145 


145 


1.7 


124 


120 


118 


117 


1.8 


95 


92 


89 


87 


1.9 


69 


67 


65 


62 


2.0 


48 


47 


45 


43 



TABLE I: Original and FSI-corrected widths in MeV for the 
Trbi mode of the tti decay, for various values of rae.x and A. 



where A is a free parameter on the order of 1 GeV. The 
strength Vq was chosen such that the squares of the po- 
tentials in Eqs. H31|) and (|34|l were equal, when integrated 
over p and q. We see that the results obtained using this 
toy potential with a comparable strength are quite simi- 
lar to those for our full model. 



rriex [GeV] 


Trbi 




1.4 


83(85) 


9(9) 


1.5 


145(149) 


7(7) 


1.6 


142(144) 


5(5) 


1.7 


118(117) 


3(4) 


1.8 


91(89) 


2(2) 


1.9 


66(64) 


1(1) 


2.0 


46(44) 


1(1) 



TABLE III: FSI-corrected widths in MeV for the yrfei and vrp 
modes of the tti decay, for various values of mex and for A = 1 
GeV (2 GeV). 




FIG. 10: FSI-corrected width of tti nbi in MeV as a 
function of the potential strength /, for rUex ~ 1.6 GeV and 
A = 2 GeV. 



mex [GeV] 


no FSI 


A = 0.5 


A = 1 


A = 2 [GeV] 


1.4 


3 


8 


8 


8 


1.5 


3 


6 


7 


7 


1.6 


3 


4 


5 


6 


1.7 


2 


2 


3 


4 


1.8 


2 


1 


2 


3 


1.9 


2 


1 


1 


2 


2.0 


1 


1 


1 


1 



TABLE II: Original and FSI-corrected widths in MeV for the 

TTp mode of the tti decay, for various values of niex and A. FIG. 11: FSI-corrected width of tti TTp in MeV as a func- 

tion of the potential strength /, for niex = 1.6 GeV and A = 2 

T„ r-n_Ul_ ITTTI J. _ :Jj.U_ £ 1 J.U 3 GeV. 
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IV. SUMMARY 

We have found that only unnaturaUy large potentials 
in the final state may change the widths for exotic meson 
decays calculated in microscopic models. In our model, 
the partial widths tend to change by only a few MeV. 
However, the irp channel is subjected to a rather large 
relative correction. It may be caused by a large value of 
the original width for the irbi mode. The net width for 
the TTi — > Tip channel is still of order of few MeV which 



agrees with the results of Refs. [11 or [13, but disagrees 
with results of Ref . . 
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